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ABSTRACT 



OPTIMIZING HULL STRUCTURAL DESIGN 
WITH 
ALUMINUM AND HYPOTHETICAL MATERIALS 

by 

HOWARD LAWRENCE SIPPLE, JR. 

LCDR. , U.S.N, 

NAVAL ENGINEER AND S,M, THESIS 

M.I.T. , 1970 

This thesis performs a modification of a computer pro- 
gram for the design of a midship section. The modification 
allows materials other than steels to "be used. The previous 
version of the progra.m had been extensively tested and proven, 
using both medium steel and high strength steels. To exer- 
cise and test the modified version of the program, various 
hypothetical combinations of modulus of elasticity and yield 
strength were used as input "materials," Test runs ware 
also made with the properties of aluminum, Tests were con- 
ducted with steels to confirm that the previous results 
could be duplicated after the modifications had been made, 

The program was formulated under merchant ship classing 
criteria of the American Bureau of Shipping, The modified 
version reflects changes, current to 1969 ■ made to the 
A.3.S, Rules since the last version of the program was writ- 
ten. The logic and limits of the program should all^w sub- 
stitution of Naval standards or other systems of classing. 
This author does not agree that this can be "easily done", 
as was previously stated, 

The program is large; it will not fit the VI AT PC R com- 
piler. The source listing for the program is retained, by 
the Department of Naval Architecture and Farine Engineering 
at M.I/T. 

Documentation of the program, detailed, instructions 
for preparation of input data, sample data and. sample com- 
puter output are included as appendices to this thesis, 
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NOMENCLATURE 

All symbols are defined locally. For the convenience 
of the reader, this nomenclature is collected below in 
alphabetical order. Greek symbols used are listed separate- 
ly, in the order of the first letter of their Engish spel- 
ling. All symbols may have slight local variations of mean- 
ing or local units. 

a Span of plating, usually in direction of applied 
loading (in) 

b Span of plating, usually perpendicular to direc- 
tion of applied loading (in) 

c A constant or a boundary condition coefficient 

c & Aluminum costs 

c Steel costs 

C A constant 

C mo Modified "old" cost ratio 

C n "New" cost ratio 

C Q "Old" cost ratio 

E Modulus of elasticity (psi) 

E' Prime symbols have been used throughout to signify 
new materials, other than steel, 

f a Fabrication ratio; 10^ greater than steel rate 

f Fabrication ratio; rate for steel 

F^ Factor which represents ratio of c^ (allowable 
stress) to a 

F T Flange thickness of I beam (in) 
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NOME NCLATURE (CONTINUED ) 

F w Flange width of I beam (in) 

Gj_ Factor which represents ratio of ~c* (allowable 
stress) to c u it 

h Water head (ft) 

I Moment of inertia, (usually in 2 ft 2 ) 

k' A boundary condition coefficient 

K Edge support condition; K = 1 for simple sup- 
ports, K = h for clamped supports 

KEY Print control symbol; KEY = Print all, KEY = 1 
Print summary only 

M Bending moment (ft-lbs or in-lbs) 

M.I.T. Massachusetts Institute of Technology 

N A.B.S. plate factor, N = cbh/12 

psi Pounds per square inch — — 

s Span (ft) 

t Thickness, usually of plates (in) 

t D Thickness of Mariner deck plating (in) 

t s Thickness of Mariner bottom shell plating (in) 

w A distributed load (lbs/ft or tons/ft) 

W n "New" material weight ratio 

W* "Old" material weight ratio 



W D Web depth of I beam (in) 

Wm Web thickness of I beam (in) 

Z Section modulus (in ft or in^) 
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NOMENCLATURE (CONTINUED ) 

Greek symbols used: 

a Used as a proportionality symbol 

6 Deflection (usually inches) 

Y a Weight of light materials in "semi -aluminum" series 

Y n Weight of "new" material 

Y Q Weight of "old" material 

Y S Weight of heavy materials in "semi-aluminum" series 

X Poisson ratio 

it Standard meaning; «n = 3.1^159 

c A stress (psi) 

c Primary stress, ship bending stress (psi) 

cr ? Secondary stress, plate-stif fener, girder bending 
stress (psi) 

00 Tertiary stress, plate bending stress (psi) 

o-.-do Apparent allowable stress implied by A.B.S, Rules; 
= 20,200 (psi) for medium steels 

Ga£ S (x) Allowable stress, comparable to o^bs ^ or steel, 
for any other material 

c_„ Critical stress for buckling (psi) 

g?4 Allowable stresses for design (i = 1, 2 or 3) (psi) 

Orn-H Ultimate strength (psi) 

c y Yield strength (psi) 
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CHAPTER 1 
INTRODUCTION 

Throughout man's history, whenever he has created a need 
to transport large loads, he has resorted to a water borne 
vehicle. Unable, himself, to carry much more than 50 - 100 
pounds for a distance of about 20 miles per day; man has con- 
stantly searched for easier ways. Land transport vehicles 
have progressed from mule back and horses through carts and 
wagons to todays' modern trucks and trains. Cargo hauling 
capabilities have increased from less than a ton carried 
about 30 miles per day to 15 or 20 tons per truck or 20 - kO 
tons per rail car over distances of about ^00 miles per day. 
Coupling 20 or 30 mules into a train or several hundred 
modern rail cars into a train multiplies the total load car- 
ried, of course; but the totals still fall far short of those 
possible with water transportation. In 1550 B.C. stone 
oblisks weighing 700 tons were transported on the Nile, [1]* 
This tonnage far over-shadowed all contemporary land trans- 
port of the day; as the almost unbelievable tonnages of the 
modern super tankers and giant bulk carriers overshadow their 
contemporary competition. 



* Numbers in square brackets refer to References at end of 
paper. 
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Ship ox-mer's demands have always pushed the naval 
architect toward greater load carrying capacity. One of the 
many decisions affecting the success of every ship design is - 
What percentage of the total ship weight must be devoted to 
self-supporting and load-supporting structure? 

The ever present demand to reduce the percentage of to- 
tal weight devoted to structural members has recently led to 
increased use of both high strength and low weight materials 
in ship design and construction. High initial cost of the 
special purpose steels and of the light weight metals, such 
as aluminum, have, however, impeded their wide spread use. 

The earliest and most extensive uses of these materials 
have occurred in "special use" vessels; where the cost con- 
siderations have been over-ridden by the uniquiness of the 
tasks that the vessel is designed to fulfill. Military ships 
fall in this category and the United States Navy has been 
among the leaders in utilizing light weight and/or high 
strength materials, especially for superstructures. The 
phenominal grox*th of complex electronic equipment, with the 
associated antennas which must be located high in the ship, 
has been one of the main factors driving designers toward 
lighter superstructure and topside weights. 

In merchant vessels, the initial use of light weight 
materials, aluminum being the one of most widespread use, 
was limited to outfit and furnishings. Eventually ladders, 
exteriors of stacks, masts, and entire superstructures were 
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constructed of aluminum, The list below indicates estimates 
of weight savings possible in various areas. The saving is 
expressed as a percentage of the equivalent steel weight. [2] 

Superstructure 35^ Q f steel weight 

Davits and life boats l±7% 

Outer stack 14.8% 

Gratings 4.9^ 

Heating coils \\yfa 

Piping, ladders, hatch covers, 31 - 36$ 
windows 

With increasing developments in metallurgy, the use of 
high strength materials in hulls, as well as superstructures, 
has been actively explored and utilized. In the case of high 
strength steels, and most other metals, the specific weight 
of the alloy is almost unchanged from that of the base metal, 
when strengthening by alloying. The weight savings are real- 
izable due to reductions in the sizes and thicknesses of 
shapes and plates needed to support the specified loadings. 
Unfortunately most applications, to date, have been hampered 
by the various existing classification societies' rules. 
These rules, in general, apply to conventional construction 
of conventional hull forms with conventional metals, and they 
do not permit full realization of the new materials' poten- 
tials. 

Any attempt to date the interest in such lightening of 

structrual weight would be pointless. One well documented 

paper presented in 19^2 to the Society of Naval Architects 
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and Marine Engineers dates the level of professional interest 
at that time, [3] Since that time many papers have been 
published on the subject. (References [2, 3, 4, 5 and 61 
represent only a very' small percentage, but were- of particu- 
lar interest to the author) 

In the field of entire hulls of light weight metals, 
specifically aluminum, a few specific dates are worthy of 
mention. In 1894 Yarrow and Company, an English firm, built 
"La Foudre," a 60-foot aluminum torpedo boat, for the French 
Government, Despite the high cost, this vessel was built 
entirely of aluminum. Thicknesses were fifty percent in ex- 
cess of those used in similar boats of steel construction. 
Two and one half tons of weight were saved; in a vessel that 
weighed only ten tons when complete, with machinery. This 
was a savings of 2$% and contributed greatly to the vessels' 
attainment of 20,5 knots average speed during full-power 
trials. This represented about a three and one-half knot 
speed advantage over contemporary craft of the same type. 
The hull was not able to withstand the marine environment, 
however, since an aluminum-copper alloy had been used. [7] 
This poor corrosion resistance and the high cost were to pre- 
vent aluminum from achieving wide acceptance in the marine 
field for many years. 

The United States Navy experimented with limited use of 
aluminum in two torpedo boats in 1899 • Great use was not 
made of this light metal until 1926, however. Under the 
Washington Disarmament Treaty Limitations the displacements 
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of capital ships was restricted. Hence the weight savings 
possible with aluminum assummed new importance. New con- 
struction and major conversions from that time to the pre- 
sent have made extensive use of aluminum, except during the 
war years when aluminum usage was restricted. Usage of 
aluminum for hulls of any appreciable size did not appear un- 
til roughly the 1950' s. During 19^8 - 1950 four aluminum 
hulled PT boats were built for the U.S. Navy. These ranged 
from 89 to 105 feet in length, and some were of riveted, 
welded and mixed riveted-welded hull construction. The Navy 
is presently receiving from contractors all aluminum motor 
boats (PGM class) of 16^ feet in length. Many smaller 
vehicles such as landing craft and other amhibious craft are 
also being built of aluminum. 

Special purpose craft being built for the Navy, includ- 
ing ground effect machines and skimmers and experimental 
hydrofoil craft, make extensive use of aluminum. The larg- 
est hydrocraft in use by the Navy is 220 feet in length, ^0 
feet in beam and displaces 300 tons. Over 110 tons of 
aluminum were used In construction of this vessel. 

Modern destroyer-type ships are being designed with ex- 
tensive topside aluminum components that result in about 150 
tons of weight savings over comparable steel designs, Many 
designers and naval architects are advocating an all-aluminum 
destroyer, predicting a stronger, lighter and faster vessel. 
The prediction is based on taking half of the expected weight 
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savings and putting it back into the ship as a stronger hull, 
still lighter than the steel hull and thus capable of greater 
speeds with the same power plant. [8] 

In the merchant marine the same trend toward larger and 
larger all-aluminum hulls is evident. The largest now in ser- 
vice is believed to be the "Sacal-Borincano. " The ship was 
placed in service in June of 19^7 . It is called - "the larg- 
est all -aluminum ship ever built" and - "the first large 
ocean-going commercial vessel with an aluminum hull," Bene- 
fits claimed over a steel design include - "reduced draft by- 
two feet, increased speed by one knot, savings in fuel, 
improved stability, and maintenance dramatically reduced. "[9] 

Thus ALUMINUM seems to have come of age as a ship build- 
ing material; and a need exists for a midship section design 
program that will handle materials other than steels. 
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. CHAPTER II 
HISTORY OF THE MIDSHIP SECTION DESIGN PROGRAM 

The history of any computer program starts with the 
pre-establishraent of a manual system for accomplishing a 
particular task. There have been many ways of deciding the 
scantlings of the midship section - the earliest was presum- 
ably a good guess which proved capable of sailing the seas 
and returning. It should be emphasized that the "returning* 1 
was of primary concern and this requirement rapidly developed 
into "return with a profitable cargo," Hence, we define 
again the designer's controversy: the hull structure must 
have enough strength (which is also to say weight) to sup- 
port the ship, with cargo, etc., at sea, in all weather, and 
yet must allow as much reserved weight and space as possible 
for cargo. 

Over the years of trial and error, a considerable his- 
tory, both of successes and failures, has accumulated. From 
this accumulation certain thumb rules can be extracted as 
being vitally necessary to a successful design. Review, 
evaluation and codification of such data is the work of the 
various classification societies. "Rules for Building and 
Classing Steel Vessels" by the American Bureau of Shipping 
(called A.B.S. Rules) is one example of such a codification, 
[10] This comprehensive body of data has for many years 
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formed the lower limit measure of structural adequacy based 
on actual operational experience. 

Unfortunately, as occurs in almost all instances when 
rules are published, there are a certain number of cases that 
fall outside the range of tabulated experience; and outside 
the range to which the rules may sensibly be extrapolated. 
As newer materials and construction techniques are tested and 
found suitable for structural use, the need for wider and 
more comprehensive design techniques is evident. This is not 
to say that the wealth of "past experience" should be dis- 
carded or ignored; but rather it should be tempered with the 
knowledge of nex-jer materials and their capabilities and 
extended on the basis of elementary and analytic concepts. 

The history of the attempt, at Massachusetts Institute 
of Technology (M.I.T.), to establish a design procedure 
based on such concepts is concerned mainly with the work of 
Professor J. Harvey Evans, of the Department of Naval Archi- 
tecture and Marine Engineering. 

In 1958 Evans presented a paper to the Society of 
Naval Architects and Marine Engineers (SNAME) entitled - 
"A Structural Analysis and Design Integration. . .with 
Application to the Midship-Section Characteristics of Trans- 
versely Framed Ships." [11] This paper was an attempt to 
establish a rational design method for the longitudinally 
effective structural material of a ship's hull. M.L. Sellers, 
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in his discussion of Reference [11] , summarized the feelings 
of many Naval Architects and others in the marine design and 
construction field. The first, often quoted, paragraph and 
the last paragraph of Sellers' discussion are quoted below: 

"When introducing young designers to the art 
of ship structural design, in order to make them 
have the proper respect for their work, we usually 
go to great lengths to impress upon them the mag- 
nitude of their task "by pointing out that a ship 
is an elastic structure weighing thousands of tons 
which is driven at a speed of many miles per hour 
through ocean waves. As such, the structure is 
subject to large static and dynamic forces and 
operates in a corrosive medium and has to last for 
many years which means, in short, that a modern 
ship is an engineering marvel. After this build- 
up, we know there is a tremendous let-down in the 
minds of these fellows when we produce the rules of 
the classification societies and show them that 
scantlings for these engineering marvels can be 
determined about as easily as reading a time table, 
Tha,t such a situation does exist is of course a 
tribute to the great work done by these societies," 

"In this paper the author has given us another' 
tool to assist in bridging the gaps between theory, 
experiment, and practice. While there appear to be 
some inconsistencies in the rules, the paper does 
demonstrate that scantlings which are satisfactory 
for service can be expressed in the terms of simple 
formulas, Once this is done we have a basis for 
further analysis and extrapotation, " 

Certainly not all discussions of the paper, in fact not 
all of Sellers' discussion, were as complimentary to the 
author. Each discussion did, however, in one form or another, 
express much of the sentiment of Sellers' last paragraph, 

"Once this is done we have a basis for further analysis 
and extrapolation." This was the vital step necessary for 
computerization also. Once the intricacies of the A.B.S. 
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Rules were reduced to formulas of proven match to the rules, 
the stage was set for eventual computerization of the 
lengthy, trial and error, iterative process of finding "the 
best design." 

Briefly the process suggested is one of determining, 
through instability criteria and stress formulas, the minimum 
thickness of hull plating and other effective load carrying 
longitudinal structure. Once an initial sizing of each 
effective member is determined, a section modulus may be com- 
puted and checked for adequacy. Modifications may then be 
made to appropriate scantlings to attain a satifactory sec- 
tion modulus from the viewpoint of hull girder bending. 
Many, many cycles or iterations are necessary to find the 
optimum solution. These iterations are the exact endeavor 
for which a computer solution is ideal. 

From Evan's 1958 paper have sprung many studies in the 

field of systematic analysis of ships' structures. At M.I.T., 

under Evans' direction, such studies have resulted in a 

detailed design computer program that has been well tested 

and proven. Many persons have contributed to this effort 

and only a few of them are represented in a listing of the 

major steps of the programs' evolution. To these unsung 

contributors is due a fitting measure of credit. The major 

milestones of the program itself are these: 

1) Weight-Strength Analysis of Cargo Ships' Structures 
by Dan Khoushy, I962 . [12] 
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. 2) Optimized Design of Midship Section Structure 
by J. Harvey Evans and Dan Khoushy, 19 63. [13] 

3) Computer Applications to the Design of a Combination 
Framed Midship Section, by Charles Roth and Donald 
Liu, 1966. [l4] 

k) Optimization of Primary Hull Structure with Mixed 
Framing Systems, by Robert D, Rockwell and Otto 
P.J. Jons, 1967. [15] 

5) Optimizing Hull Structural Design with High Strength 
Steels, by J. Harvey Evans, I967. [16] 

Each of the above have modified and updated and/or test- 
ed the then current version of the computer program. Ini- 
tially designed to handle only transverse framing systems, 
the program will now handle transverse , longitudinal and 
mixed framing systems. It has been used in exhaustive 
parametric studies to define optimum least weight and least 
cost designs and has been used to determine the optimum 
distribution of various combinations of high strength steels. 
[12, 13, 14, 15 and 16] Not until this work, however, have 
materials other than steels and materials with varying 
elastic modulus been investigated with the program. 
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CHAPTER III 
MODIFIC ATIONS TO THE COMP UTER PROGRAM 

General Comments on Modifications 

As stated In Chapter II many authors' variations have 
contributed to the present version of the program. Many- 
statements in previous theses, concerned with the midship 
section program, have indicated that the program was written 
in very general terms. As initially envisioned, any modifi- 
cations to the program were to be a very minor part of this 
thesis. The main effort was to be directed toward an econo- 
mic study of the .cost effectiveness of utilizing aluminum as 
a hull material, knowing that a significant premium must be 
paid for material cost. The question to be answered was - 
Could the optimum midship section design save enough hull 
structural weight to pay off the material cost premium in a 
reasonable length of time and then earn the owner a larger 
profit? 

It was found that the program was functionally designed 
to be general but that realistically, as written, it could 
handle only steels. Many instances were found where vari- 
ables such as specific weight were written in as constants, 
(0.283 pounds per cubic inch for steels) rather than as a 
variable from the array of material properties. There is no 
record whether this was an expedient to meet some past thesis 
deadline or merely fore-knowledge that only steels would be 
used in the immediate future, 
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The section that designs the transverse structure, which 
was a separately designed addition to the program, was based 
entirely on medium steel formulas, with pre-computed con- 
stants determined by fitting the formula to the American 
Bureau of Shipping (A.B.S.) Rules. These formulas were 
constructed in a manner that made no allowance for variations 
In properties, modulus of elasticity or yield strength, from 
the values for medium steel. 

Although the program was written to handle the design of 
the midship section using more than one metal and although 
the program has been successfully utilized in this fashion, 
the success achieved depended on the fact that all the com- 
binations of materials used were steels. (Reference [l6] 
used three sections Bottom, Side and Top, or main deck, with 
locations of the interfaces a variable) Combinations of 
other materials, or of steels with other materials would not 
have worked. 

Early in the preparation of this thesis, it became ap- 
parent that the modification of the 75 subroutines of this 
program, amounting to almost 6000 Fortran cards, was a task 
of signicantly larger proportions than had been previously 
anticipated. As a result the portion of the time budget 
originally scheduled for the cost effectiveness analysis 
of the results of this study was re-directed toward the 
thorough modification of the program to enable any material 
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or combination of materials to be utilized. Another tangent 
was added at this time also. It was decided to run a para- 
metric study of the effect of variation of combinations of 
modulus of elasticity and yield strength. 

The "hypothetical materials" that were "constructed" 
have the following properties - modulus of elasticity varies 
from 10 x 10 6 (psi) through *K) x 10° (psi), in steps of 
10 x 10° (psi); yield strength varies from 30 x 10^ (psi) 
through 150 x 10^ (psi), in steps of 40 x 10- 5 (psi). These 
various hypothetical materials present a matrix of sixteen 
possible combinations; for each of which the "optimum" frame 
spacing is desired. 

One of the major obstacles encountered in attempting to 
modify the old version of the program was the myriad of 
cryptic abbreviations used as Fortran names for the variables. 
In order to modify the program Intelligently it was first 
necessary to dissect it, in the smallest detail, to learn 
not only what each variable name meant but also what each 
element of each variable arra3' contained. This effort has 
caused the author to recommend that a tabulation of the sort 
found in Appendix 2 be a requirement for future computer 
programs; especially those that are liable to be expanded 
and modified. This subject is discussed in more detail in 
the section on recommendations for further study. 

The remainder of this chapter is divided into two parts- 
Modifications to Data Deck Setup and Modifications to Accomo- 
date Variable Material Properties. 
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Modifications to Data Deck Setup 

The previous program (program version 5) is Ideally 
setup for parametric studies. The data input system is 
extremely flexible allowing almost all variables to be 
varied at will. Any of the data input sections ($MS ******) 
can be changed individually. (With due regard to keeping the 
number of materials in the $MS PROPERTIES OF MATERIALS 
section equal to the number of cost sets in $>MS COST 
section and having $MS COST data always follow the 
$MS PROPERTIES OF MATERIALS data in the data deck setup) 
Only two major changes were made in the data deck setup; and 
one of these was more on the order of an addition to the 
system, rather than a modification. 

A separate data section "PROPERTIES OF MATERIALS" was 
created and a new subroutine, called PROMAT, was constructed 
to read in the properties of the various materials that the 
program would utilize, The old "MATERIAL" section was renam- 
ed "MATERIAL ASSIGNMENT" and presently handles only that 
explicit task; reading in the assignment of particular 
materials to the various sections of the ship. 

The other major change concerns an Indicator to control 
the printed computer output. Basically this indicator per- 
mits either all detailed Intermediate results and the summary 
or only the summary of results to be printed. 
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A separate section, Appendix 3, explains, in step-by-step 
fashion, the construction of a sample data deck. Appendix k 
contains .a copy of a sample data deck. 

The reasons and mechanics of the above changes, and a 
few other minor changes, follow: 

In the previous program the properties of the metals to 
be used during computation were inserted in the Fortran 
source deck of the program as DATA statements. They could be 
changed at will by simply exchanging new cards for old. The 
length of the program and the resultant cost for reading in 
and compiling the Fortran source decks for each data run, 
led to the conclusion that operation with a compiled binary, 
or object, deck would be much more efficient. This certainly 
proved to be the case. 

To permit the changing of the materials while operating 
with the binary deck, a new subroutine (PROMAT) was construct- 
ed to read in the properties of each material desired and 
from these properties to compute the allowable stresses to be 
used later as design limits. The subroutine data input for- 
mat for PROMAT is detailed in Appendix 3. 

The subroutine PROMAT expects to read first the number 
of materials to be read in. This number is limited to five 
or below. This same number is used in the INCOST subroutine 
for reading in cost data for each material. Because of this 
the material data (#MS PROPERTIES OF MATERIALS) must always 
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precede the cost data ($MS COST DATA) In the data deck 
setup. The number of sets of cost data (cost of plates, cost 
of stiffeners, and labor factor) submitted must alxrays match 
the number of materials submitted as data. 

PROMAT next expects to read two cards for each material. 
The first card contains the material type number, for identi- 
fication, the yield strength, the ultimate strength and an 
allowable stress factor. The second card contains the modu- 
lus of elasticity, Poisson ratio and the weight of the mate- 
rial. The allowable stress factor and Poisson ratio are 
diraensionless. The weight must be in units of pounds per 
cubic inch. All other properties are in units of pounds per 
square inch. 

The allowable stress factor permits two methods of se- 
lecting the allowable primary, secondary and tertiary stress- 
es. Allowable stresses are called a-, f a ? and 031 as defined 
in [11] t with the usual meanings: 0, = ship bending stress, 
a? = girder bending stress (plate-stlffener combination), 
and Crt ~ plate bending stress. 

The previously published and generally accepted figures 
of cr = 20,000 (psi), c 2 = 27,000 (psl) and = 32,000 (psi) 
for a medium steel, with a yield strength of 32,000 (psi) 
and an ultimate strength of 60,000 (psl), were used as a 
basic starting point. The numbers, listed as "Based on a 
and cr U ] t " in Table 3 of Reference [16] , were obtained by 
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Evans in the following manner, A factor relating each sigma 
stress to the yield strength and another factor relating 
each to the ultimate strength were defined: 

0i cri 

F, = G 



i w l " 

c y G ult 



F-, is therefore defined as: 



0-, 20000 

F-, = _ = = 0.625 

c 32000 



The values of F and G are: 

F^^ = 0.625 G = 0.3333 

F 2 = 0.84-375 G 2 = 0.4-50 

F3 = 1.0 G 3 = 0.5333 



Using these values, a based on c and o. "based on 
Cy.-]*. can be found: 

*i = Vy \ = G i0 u lt 

Thus two limiting values of allowable stress are gene- 
rated; one based on <j v and one based on cr n .. There is lit- 
tie agreement as to which should be used as the base for 
limiting design stresses, 

R.G. Kline has dealt with high strength steels, under 
brittle fracture and fatigue conditions. [k~] He concluded 
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that acceptable allocable stresses for high strength steels, 
with a y = 50,000 (psl) and 100 ,000 (psi), could be generated 
by simple .multiplication of the accepted allowable stresses 
for medium steel by 1,3 a nd 1,8, respectively. 

The allowable stresses generated by the above methods 
are shown in Table 1, This table is a reproduction of Table 
3 of Reference [16], To it have been added the various alter- 
natives which may be generated by the present (version 6) 
program. 

In general, the allowable stresses generated by Kline 
fall roughly mid-way between those based on c and c ult 
The values used for calculations in Reference [16] were 
Kline* s values with ?;ome small changes made in the allowable 
tertiary, or plate stress, values. It was decided to con- 
struct the program in such a manner that the user could mod- 
ify the medium steel allowable stresses by either a simple 
multiplication factor, such as Kline proposed, or have the 
program generate the average between the ratioed values 
based on c and a u ^^.. 

The "Allowable ;;tress Factor" listed as an input in the 
$MS PROPERTIES OF MATERIALS data section is used to select 
which method will bo used to set allowable stresses within 
the program. If thu Allowable Stress Factor is set equal to 
1.0 then each medium steel allowable stress is multiplied by 
1.0 and the medium nteel allowable stresses result. If any 
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other positive (non-zero) Factor is set in as input, the 
same sort of multiplication of the medium steel allowable 
stresses by that Factor occurs, to generate new allowable 
stresses for the new material. If the Allowable Stress 
Factor is set equal to zero, the program utilizes the pre- 
viously stated values of F, and G. and the c and r> , . values 

11 y ult 

of the material to generate the average of o^ based on q 
and c^ based on c^f These averages then become the allow- 
able primary, secondary or tertiary stresses (depending on 
the value of i) used in the program. 

To control the printed output, all print subroutines, 
except the summary subroutine, are set-up to be keyed by an 
input indicator, called "KEY." When KEY is set equal to 
zero (KEY = 0) all subroutines will print. This givesjthe 
user all the detailed intermediate results as well as the 
summary of results. When KEY is set equal to one (KEY = 1) 
only the summary subroutines will print. This summarization 
always includes the print out of final plate thicknesses. 
Plate thicknesses were considered desirable results even if 
only the summary printout is requested. 

When parametric studies are being conducted the results 
desired are provided by the summary. Total weight and cost 
per foot, and weight breakdown into transverse and longitudi- 
nal structure are printed. Longitudinal weights are further 
broken into weights of various materials when more than one 
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material is used. Cost is broken into material cost and 
labor cost subtotals, 

Two ratios are provided also by the summary. Thes^ are 
a weight ratio and a cost rati.o and are formed with input 
base weight and base cost. Most studies done at M.I. T. to 
date have used the Mariner, hull j and the base weight and cost 
figures used represent this hull: transversely framed, made 
entirely from medium steel, and with a frame spacing of 2,5 
feet. 

Details provided by the remaining print subroutines 
includes coordinate locations of all plate seajns and all 
stiffeners; head loadings at each of the above locations; 
bottom structure weights: area, moment, and inertia constants 
used: initial plate thicknesses assumed; initial stiffener 
sizes assumed; final stiffener sizes; final plate thicknesses; 
transverse weight breakdown; deck and shell frame weight 
breakdown; cost breakdown. 

When designing a midship section, naturally the user 
will want all the design details printed out. For parametric 
studies of a large number of ships, or many variations of the 
same ship, however, the user will probably be satisfied with 
the summaries. This may be the case for the designer above 
also, until the optimum has been found; then a rerun "keyed" 
to print all detail would provide the "design." The differ- 
ence between KEY = 1 and KEY = is about twenty pages of 

printed output. 
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Another minor change provides for several materials to 
be displayed in the STRUCTURAL DATA input section. This 
section prints the "bending moment to be designed for, the 
limiting design stress and the resulting required section 
modulus for both bottom and upper deck areas. Previously 
only one material could be entered here. This information is 
not used within the program and is printed out only for com- 
parison with the program generated values of section modulus. 
If desired the program will generate this data using familiar 
formulas for medium steel, when no data is provided as imput. 

Another change provides for the selection of the mate- 
rial used for design of the transverse structure. In the 
previous version all transverse structure xras of medium steel. 
This made some sense when constructing of high strength 
steels for longitudinally effective material, as it would" 
probably be the practical solution to design the transverse 
structure of the cheaper medium steel. This can still be 
designed in, if desired, since the material for transverse 
structure is a separate selection made in the $MS MATERIAL 
ASSIGNMENT section. 

A data check was written into the subroutine INCOST, 
which reads in the cost data, to prevent accumulating a 
different number of sets of material properties and sets of 
cost data. This check requires these numbers to be matched 
and hopefully helps avoid costs being matched with wrong mat- 
erials or no costs being entered for some material, etc. 
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Modificat ions t o Accomodate Variab l e Material Propert ies 

The computer program has been constructed to satisfy 
the following basic design equations: 



Ch + CF 2 ^ C 2 



Sigma bar quantities ( ) represent allowable stress- 
es. Un-barred quantities ( ) represent nominal hull 
stresses, [11,13 and 16] 

These design equations are based on medium steel, where 
the ultimate strength is roughly a factor of two times the 
yield strength. Provisions were made in the modified pro- 
gram to allow the user to adjust the allowable stresses by 
ratios related both to the yield and the ultimate strengths 
of the material to be used. This modification was explained 
in the previous section since it was directly related to the 
data input of the new subroutine PROHAT. 

Further modificationr to the design equations were 
found to be necessary to permit the use of materials other 
than steels. These modifications were mainly associated 
with the section of the program that designs the transverse 
structure. In the longitudinal structure sections the design 
formulas are based on more elementary principles rather than 
directly on the A.B.S. Rules' formulas. The transverse 
structure design section is less rationally developed. It is 
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still directly related to the A.B.S. Rules and sizes trans- 
verse structural members using formulas derived from design 
loadings deduced from scantling tables of the A.B.S. , and in 
a fevr cases taken directly from the Rules. [12] 

Modific ations to the Formulas For t h e Desig n of Deck 
Beam s_ and S i de F rames 

Starting with the general bending moment equation: 

2 
ws 

M = x 12 (a) 

k 1 

M 

and the relation: o = — (t>) 

Z 

where: M = bending moment (in-lbs) 
w = distributed load (lb/ft) 
s - span (ft) 
Z = section modulus (in-^) 
o = allowable bending stress (psi) 
k' = boundary condition coefficient 
The o we are dealing with here is a 2 » the allowable 
secondary or girder bending stress, for the beam. For pur- 
poses of identification of this particular stress, in the 
work following, it will be called a A -D C » and specifically 
refers to medium steel unless otherwise noted. 

From equations (a) and (b) , Z is determined: 

2 
12ws 

k ' <?ABS 
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The A.B.S. Rules define the required section modulus 
for beams as: 

Z = 0.00315 Ns 2 ' (d) 

where: N = cbh/12 (e) 

so that: Z = 0.0002625 cbhs 2 (f) 

where: c = a boundary condition coefficient 
b = frame spacing (ft) 
h = pressure head (ft) 
Z = section modulus (in^) 
s = span (ft) 

The computer program sizes the beam under consideration 
by forming a ratio of the required section modulus, from (f) 
above, to the section modulus of a standard stiffener. This 
ratio is called the proportionality factor and tells what 
proportion of the standard stiffener, either smaller or 
larger, Is required. 

If the section modulus, Z, is inversely proportional to 
the op stress, as in (c); and if the material to be used 
has an allowable stress not equal to that of medium steel; 
then the required section modulus should reflect the differ- 
ences in allowed stresses by requiring an inversely corres- 
ponding section modulus. An allowable stress that is smaller 
than 0« would require a larger section modulus, etc. 

Let 0VBq(x) represent the allowable stress of the new 
material (x). 
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Thens Z = — ( g ») 



. 12ws 2 
3 - 


k * °ABS X 


a ABS (x) 
a ABS 


2 

12V7S 

Z - x 


a ABS 



(s) 

k ' a ABS a ABS^ x ' 

Hence the required Z for medium steel should be multi- 
plied by the ratio of medium steel allowable stress to the 
allowable stress of the new material. 

Another factor complicates this simple solution, how- 
ever. The 0Vo S defined above is not the c 2 defined in the 
PROPERTIES OF MATERIALS section, Evans and Khoushy [13] 
found that by working backwards (with the knowledge of nor- 
mal sizes for spans, construction practices, the usual range 
of stiffness and the extremes of 1.0 and 0,0 for the distri- 
bution factors for simple supports and fixed ends, respec- 
tively) the allowable stress implied by the A.B.S. Rules is 
9 (tons/in 2 ) or 20,200 (psi). 

As explained previously in the PROPERTIES OF MATERIALS 
section, this program uses the more generally accepted 
27,000 (psi) for 02 » the allowable stress for medium steel. 
As found by Evans and Khoushy: 

a ABS ( steel ) = 2 0f200 (psi) (h) 

with: a 2 (steel) = 2?,000 (psi) as before (i) 



Co (steel) 

then: cTatdq (steel) = 20,200 x — (y) 

A * 27,000 



cfo( steel) 

a ABS ( steel ) = 20,200 x z — — (J) 

<j2 (steel) 

where c 2 ^ as its previously defined meanings 

On - allowable secondary stress as computed 
in the PROPERTIES OF MATERIALS section 

Therefore, 

a A gg(x) for any other material: 

tfo(x) 

a ABS (x) = 20 » 200 x — — — (*) 

o>> ( steel) 

Now from (g) , the required ratio for modifying the sec- 
tion modulus is: 

p ABS (steel) 
a ABS^ x) 

This ratio is called Transverse NAterial RATio, or 
TMARAT, in the program. 

CTABS( steel) 
TMARAT = —— ~- (1) 

a ABS (x) 



But, from (k) : 



c 2 M 

^ABS (x) = 20 » 200 x =r — 

02 ( steel) 



(k) 



and from (h) a ABS ( steel) = 20,200 (h) 
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20,200 
Therefore: TMARAT = — — _ - (m) 

C2( x ) 
20,200 x — 

CpC steel) 



Co( steel) 
TMARAT = — (m) 

Hence, each time a beam section modulus Is computed "by 
the program it is modified by the ratio TMARAT. 
Thus (f) multiplied by TMARAT becomes; 

7 CpC steel) 
Z = 0.0002625 cbhs x _ (n) 

c 2 (x) 

Vertical side frame members are designed by the program 
in a like manner, using different formulas of course, result- 
ing in a specific required section modulus. This required 
modulus is also multiplied by the ratio TMARAT to reflect 
material property variation. 

Mo dific ations to the Formulas For Deep Web Frame Design 
A.B.S, requirements for deep web frame design are listed 
below: 

Deck Fra mes Side Frames 

W D ^.7xs ^ 1 . 5 x s 

W T ^ .01 x W D + 0.16 ^ .01 x W D + 0.14 

VU ^ .3 + .01 x F m x F TJ ^ .56 



h2 



- where: . w" D = web depth depth (In) 
Wm = web thickness (In) 
F w = flange width (in) 
F T = flange thickness (in) 
s = span of member (ft) 

0.l6 and 0.1^- additions are wastage allowances 
and were not modified 
The Rules state the requirements for web frame dimen- 
sions in terms of minimums, with no reasons given. This 
leads one to speculate that this minimum limit is a stiffness 
criterion and presumably is based on a deflection limitation. 
One such deflection criterion Is: 

s 



or: 



6 ^ — 
360 





6 = f(EI) 




ws 




6 = - — 


(for a distributed load w, 


cEI 


span s, and simple support 




conditions) 



where: 5 = deflection (usually Inches) 
s = span (ft) 

E ss modulus of elasticity (psi) 
I = moment of inertia (In 2 ft 2 ) 
w = distributed load (lbs/ft) 
c = constant 
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To require the same limitation or same maximum deflec- 
tion, when using a material other than steel, requires that: 

ws w* (s ' ) 



cEI c'E'I' 



where: prime quantities represent new material proper- 



ties. 



The end conditions, loading and span all remain constant. 



Therefore : 



EI E'l 



E 



or 



!• = I x 



The program uses an I-beam to simulate the web frame 
member. The moment of inertia of an I-beam is found: 



I about own axis = I TT ~-u + 2L lovl „ 

weu i lange 



: web = ~ W T V 



I flange ~ F t x F W x 



W 



D 



W D F T 

Assuming — — ^>> 

2 2 

w T w D 3 P T P W W D ' 
Therefore: I = + 2x L - 
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Fig. 1 
I-Beam Detail 



W* 



W T W D 3 F T F„W D 2 

r I = + 

.12 2 



W,p and F T are approximately equal 
W D and F w are of the same magnitude 

3 
Therefore: I is proportional to W Q 

I a W D 3 I' a (W D ») 3 

E 

Since: I', = I x * 

E' 

■p 

Then: ( w n' )3 a (W D^ x 

u E' 






The current program modifies W« of the A.B.S. Rules by 
multiplying by the cube root of the ratio of the steel 
elastic modulus to the new material elastic modulus. 



E 
Thus : W„ ^ . 7 x s x 3 

\J 

E' 

In the program the ratio is called EPRIE, and the 

E 

above limit on Wq becomes: 

1 

W D ^ .7 x s x - — 

\3/EPRIE 

^5 



Modifications to the Formulas For Span to Th ickness 
Ratios 

. . The program requires that the following proportions be 
maintained in web frame design: 

W D , F W 

~ < 60. — <• 10. 

W T F T 

A similar requirement is imposed on vertical floors in 
the bottom design: 

b 

_ ^ 96 
t 

where: b = span dimension of plate (in) 

t = thickness (in) 
These requirements are obviously to satisfy a buckling 
criterion. 

Ktt Et 

For buckling: a = ~ ~--^ 

cr I2(l-x 2 )b^ 

where: = critical stress for buckling (psi) 
cr 

K = constant depending on aspect ratio 

E = modulus of elasticity (psi) 

t = thickness of plate (in) 

b = span dimension of plate (in) 

X = Poisson ratio 
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Solving for 



Ktt 2 E 



CE 



t 2 l?Al-\ 2 )o CY (1-x 2 ) a cr 



b 
t 



CE 



(1-X*) c CT 



The relationship when the material is changed, where the 
prime quantities represent the new materials 



b' 
t' 



E 



(l-X^'ccr' 



E 



(1-X )a cr 



and if cr cr = c 2 



b f b 


E« c 2 (1-X 2 ) 


t* t 


E o 2 f (1-X 2 )' 


ed before : 


— - = EPRIE 



a 2 « 



= TMARAT 
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Define: 



(1-X 2 ) 

(i-x 2 )' 



= PPPRI 



Now: 



b' 
f 



(EPRIE x TMARAT X PPPRI) 



Also define: ESIGRA = (EPRIE x TMARAT x PPPRI) 2 



Then: 



b» 
f 



X ESIGRA 



Thus any such buckling ratio limit that is imposed, by 
the program is modified by the multiple ESIGRA. 
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CHAPTER IV 
ALUM I NUM,_ I NVES T IG ATIONS 

General Com ments 

The original scope of the aluminum investigation was 
drastically reduced when the magnitude of the computer pro- 
gram modification became evident. As completed, two differ- 
ent uses were made of the computer program during the 
aluminum investigations. One part attempted to determine if 
unit weight and cost predictions could be made by knowing 
the relationship of variables between a previously plotted 
material and a new material. Another part of the investiga- 
tion found an "optimum" frame spacing for a Mariner hull 
built entirely of aluminum. 

Definition of Semi -Alu minum Materials 

The first objective was accomplished by devising a 
matrix of different combinations of weights, costs of ma- 
terials and fabrication cost ratios. When making this series 
of computer runs the following properties were kept constant: 

a. All stresses were based on aluminum properties; 
C y = 19.000 (psi), cr ult = 35,000 (psi) 

b. Modulus of elasticity, 10 x 10 6 (psi) and Poisson 
ratio, 0.33» of aluminum were used. 

Table 2 shows the combinations of variables used to 
make up the eight "Semi -Aluminum" materials. 
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Table 2 Semi -Aluminum Materials' Properties 

Heavy (H) Li ght ( L) 

V Y s' V f s V V f s 

2) Y 8 . o a . f a Y a . V f a 

3) y » c » f Y« » c » f 

s a s a a s 

*0 Y » c » f Y Q » c f 
s a a a a, a 



where: Yo = weight of steel 

Y = weight of aluminum 
a 

c_ = cost of steel 

c = cost of aluminum 
a 

f = 1.0 fabrication cost ratio 

S 

f =1.1 fabrication cost ratio 
a 

The computed results of the various combinations are 
plotted in Figure 2 and tabulated in Table 12, in Appendix 1. 
The prediction value of the plot is illustrated by the fol- 
lowing: 

Suppose that point 1H represents a previously plotted 
material (1H) with the strength limitations of aluminum and 
the weight and cost of steel. A "new" material (^L) is 
found with the same strength properties but with a weight and 
cost of aluminum. 

The weight ratio of the new material to the old is 
found : 
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where: y = weight of new material (lbs/irr) 
Y Q = weight of old material (lbs/in^) 
If the weight ratio of the old material is modified by 
this factor, the predicted weight ratio of the new material 
is found: 

W n = .339 x W Q = .339 x W Q 
W = A9^ 

where: W = predicted weight ratio of new material 

W Q = weight ratio of old materia.! 
This point Is projected vertically downward from 1H and 
is labeled IP. 

Next the ratio of material costs is found. 

" Old" (steel ) "New" (aluminum) 

Plate cost .075 .532 

- ■ Stiff ener cost .075 .625 

The ratio of new material cost to old material cost for 

plates and stiff eners Is found: 

.532 

Plates - -7.09 

.075 

.625 
Stiff eners = 8.3^ 

.075 



Assuming that 75^ of the hull weight is plates and 25, £ 
is stiff eners: 
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. An adjusted cost ratio for . new to old material^ is 
found: Material cost ratio = 7.4-0 

• « In applying the cost differential ratio, it must be 
realized that the cost r^tio includes both material cost and 
labor costs. Figure 3 plots the subtotal breakdown of the 
cost data recorded for t.| 1Q "hypothetical" materials discussed 
in Chapter V. The plot shows the relationship of labor and 
material costs and the t.>tal costs. Figure 4 plots the same 
data points; each as a Percentage of the respective total 
cost. Figure 5 plots tho subtotal breakdown of the weight 
data recorded for the sr»., lG hypothetical materials. Tabula- 
tions of this data is foi m a in Table 11, Appendix 1. From 
Figure 4 it is estimated that, for a material with a yield 
strength of 19,000 (psi), the labor cost is approximately 
60% of the total. 

Labor costs in this series of studies has been inflated 
by 100^ to reflect overh. ?ad . ln other words tv/ice the average 
salary has been used to compute the labor cost. Therefore 
the true labor cost is ap proX imately one third of the total 
cost, overhead is one third and materials make up the remain- 
ing third. In the discu-. sionf however, the two thirds which 
comprise the labor costs and the overhead will collectively 
be called - labor costs. 

Using the labor costs as 60% of the total, and consider- 
ing these costs as constant the relation of the predicted 
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new cost to the old cost is as follows: 

Labor costs = .60 x C (Labor costs remain constant) 
Material costs = .4-0 x C Q (is modified by cost ratio) 
Therefore: Total cost (or total cost ratio) - 

Labor costs + Material costs (modified) 

C n = C Q (.6 + A x 7.^0) 
C n = 1.158(3.56) 

where: C n = cost ratio of new material 
C Q » cost ratio of old material 

This predicts the cost ratio of point 3H with fair 

accuracy, however the desired prediction is for a lighter 

material. To correct this it must be recalled that both 

materials have the same allowable stresses, hence the volume 

of structural materia.1 will be the same for each. Since the 

weights are different, however, the total weight will be 

different. Cost figures are in terms of $ per pound, so this 

difference in weight is reflected in the cost ratio. To 

reflect this change, k-0% of the old material cost ratio, 

representing the material cost, should be reduced by the 

weight ratio of the materials. The modified old cost ratio 

is: 

C mo = .k x 1.158 x .339 + .6 x 1.158 

C mo = -^7 + - 6 ?5 
c mo " « 8 ^ 

where: C = modified cost ratio for old material 

mo 
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Now using the same prediction formula as before, a cost 
ratio for the new light material can "be found: 

C n = C mo ( ^ x 7.^0 + - 6 > 
c n = c mo (3 --5 6) = - 8 52(3.56) 
c n = 3.03 

This translates point IP, horizontally to point 2P. 
This cost ratio Is considerably higher than the actual light 
weight cost ratio (4-L). By this time many errors have accum- 
ulated due to the various assumptions and approximations 
made. Never-the-iess the predicted point does show the 
order of magnitude of the weight saving achieved and the 
cost premium paid. Refinement of the above estimates by 
further experimentation should produce much closer predic- 
tions, 

Predictions for another material are compared with a 
steel performance norm curve in Figure 1*4-, in Chapter V. 
Such comparisons could predict whether the new material 
would be competitive with steels as a ship construction 
material. 

M3^II^L-£bM^]^ _. Frame Sp acing and Aluminum vs. S teel Econo- 
mic Analysis 

The second objective, that of determining the optimum 
frame spacing for an all aluminum hull, was accomplished by 
making several runs at various frame spacings. The results 
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are. plotted in Figure 6. and tabulated in Table 1^ in Appen- 
dix 1. The optimum frame spacing is 28 inches. 

. From this optimum an estimate can be made of the total 
weight and cost of the aluminum hull. The difference in 
total weight between a steel hull and an aluminum hull 
defines the extra revenue cargo that could be carried by the 
aluminum ship due to hull weight savings. The total cost 
difference defines the extra premium that must be paid for 
the higher cost of aluminum. To be economically attractive 
the aluminum hull must be able to carry enough extra revenue 
cargo to payoff the initial cost premium within a reasonable 
time, and then earn a higher profit for her owner for the 
remainder of the ships* lifo cycle. 

Evans has reported - "n reasonably consistent relation- 
ship (exists) between this unit weight and the total hull 
steel weight (expressed as n 'prismatic weight coefficient'), 
at least within a particular ship type." [16] For the pur- 
pose of this brief economic analysis it will be assumed that 
such a relationship exists, both for steel and for aluminum. 
It will be further assumed that the same relationship exists 
between the unit cost and the total cost. The relationship 
assumed will be simply that the unit weight, or cost, multi- 
plied by the vessel length will produce the total weight, or 
total cost, respectively. Xany arguments can be immediately 
advanced against each of these broad assumptions. The author 
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is satisfied with only the first; a relationship between 
unit weight and total hull weight, and even here it is not 
suggested that the relationship is actually as simple as the 
one assumed. Such simple assumptions will serve however to 
Illustrate the point. 

Using the above assumptions, the following figures are 
found : 

Mariner hull length L = 528.5 (ft) 

Steel 
Steel Aluminum Ratio — 



Aluminum 



Total Hull Weight (tons) 3295.73 1622.50 2.03 

Total coat ($ x 10~ 6 ) 1.4208 2.8358 0.501 



Weight Difference 1673.23 (tens) 

Weight Savings 50.8^ of steel hull weight 

Cost Premium $1.4-150 x 10 6 

For the Mariner class vessel, the following yearly 
cycle is assumed: 

1) One round trip (assumed 3000 mile voyage) per 
month for a total of 12 round trips per year, 
and 24 cargo trips. 

2 ) Carry cargo both ways ; in order to take full 
advantage of weight saved, consider fully loaded 
on all trips. This assumption is very bad for 
this kind of cargo ship but makes good sense for 
bulk carrier types (one way only) 

3) Cargo Freight rate of $30 per ton. 
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Increased Revenue Per Year (IRPY) = Number of trips x 

Increased tonnage per trip x Freight rate per ton 

IRPY = Zh x 1673.23 x 30 
= $1.20^7 x 10 6 

Thus, with the very questionable assumptions made, the 

total aluminum construction premium could be recovered within 
1.09 years. 

This very naive economic analysis is based on the fol- 
lowing properties for aluminum: 

Yield strength (psi) 19,000 

Ultimate strength (psi) 35,000 

Modulus of elasticity (psi) 10 x 10 6 

Poisson ratio 0.33 

Weight (pounds per cubic inch) 0.096 

Fabrication cost ratio 1.0 

(same fabrication costs as steel) 

Cost of plates (per pound) $0,532 

Sizes: Thickness range: 0.75-0.25 inches 
Length range: 72-2^0, 2^0-360, 

360-^80,^80-612 inches 
Width range: 60-80, 80-132 inches. 

Cost of stiff eners (per pound) $0,625 

Sizes: Weighted average of circle 
sizes 6-8, 8-10 and 10-12 
inches, on a 20-^0-^0$ basis, 
respectively, 

These properties are roughly representative of the 

aluminum alloy 5086, one of the most commonly used alloys 

in aluminum marine construction. 
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The value of 19,000 (psi) used for yield strength is 
based on an average of the present data available on as 
welded strengths of both plate and stiffener shapes. This 
19 i 000 (psi) is considerably lower than the "advertised" 
yield strengths of some of the heat-treatable alloys, 
however these alloys have not yet been used extensively 
for marine structures. 

The reduction in total hull weight, as a percentage of 
the steel hull weight, agrees very closely with other studies. 
Maclntyre, in a paper concerning ore carriers [5] , reported 
that the aliiminum hull weight was $£>% of the weight of the 
steel hull; or a weight savings of kh% of the steel hull 
weight. Approximate calculations based on Maclntyres' paper 
[5] indicates that one of his aluminum ore carriers would 
repay the construction cost premium in about 12 years. The 
costs quoted in [5] are from bids" on the designs. 
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■ CHAPTER V 
HYPOTHETICAL MATERIALS I NVESTIGATION 
Definition of "Hypo thet ical Ma terials" 

As mentioned earlier one series of investigations con- 
cerned hypothetical materials. The materials were "construct- 
ed" by varying the modulus of elasticity (E) over the range 
from 10 x 10 6 (psi) to ^0 x 10° (psi) in steps of 10 x 10 6 
(psi) while yield strength (q ) varied from 30 x 10^ (psi) to 
150 x 10^ (psi) in steps of kO x 10^ (psi). One combination 
of these two properties, c y = 30 x 10^ (psi) and E = 30 x 10 
(psi) , is very close to the values for medium steels. The 
other fifteen possible combinations do not represent any 
real material. 

Throughout this paper the various "hypothetical materi- 
als" will be referred to by their yield strength value, c v 
times 10"^, over their modulus of elasticity, E times 10"" . 
Therefore the material with yield strength equal to 70 x 10-^ 
(psi) and modulus of elasticity equal to 40 x 10 (psi) will 
be called 70/^0. The name is always 

c y x 10-3 



E x 10 



-6 
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Table 3 Matrix of Hypothetical Materials (p /E ) 



C /E 



10 x 10° 20 x 10 

30 x io 3 30/10, 30/20 

70 x 10 3 70/10 70/20 

110 x 10 3 110/10 110/20 

150 X 10 3 150/10 150/20 
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30 x 10 6 kO x 10 6 



30/30 
70/30 
110/30 
150/30 



30A0 

70/40 

110A0 

150A0 



These hypothetical combinations of modulus of elastici- 
ty and yield strength were investigated with two main objec- 
tives. First, and. most logical, they would exercise the 
modified program in handling variations of E and c v . Second, 
It was anticipated that, once an optimum frame spacing was 
obtained for each trial "material," a grid could be plotted 
which would enable a prediction of the optimum spacing to be 
made. More important, since a previous study had delineated 
an attainable limit curve, or performance norm, for steels, 
(least weight, least cost studies for various steel combina- 
tions: Reference [l6])the same grid could predict whether 
the new material could be competitive with the steels. If 
certain frame spacings resulted in plots below the steel 
limit curve then the new material could be considered com- 
petitive with steels. See Figure iM, 

All plots are a non-dimensionalized, weight per foot 
ratio versus a non-dimensionalized, cost per foot ratio. 
The basis for the ratio is, in all cases, the weight per foot 



65 



and cost per foot for the Mariner all medium steel hull, 
transversely framed, and with a frame spacing of 2,5 feet. 
Subsequent results are divided by these base values (6.236 
tons per foot and 268^.7 dollars per foot) to form the ratio 
for plotting. 

Determination of Trial Frame Spacing s 

In order to plan the frame spacings to be used in the 
search for the optimum spacing, several attempts were made to 
estimate optimum spacing. 

First from Sezawa: 
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the following table can be generated. 
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Table 4 
Optimum Frame Spacing to Plate Thickness R atio (a/t ) 

C /E 10 xlO 6 20 x 10 6 30 x 10 6 40 x 10 6 

30 x io3 18.25 25.85 31.65 36.50 

?0 x 10 3 11.95 16.90 20.70 23.95 

110 x io3 9.55 13.50 16.50 19.10 

150 x 10 3 8.15 11.55 14.10 16.30 

By making some reasonable assumptions about what per- 
centage various locations contribute to the section modulus, 
a relationship for t (thickness of plating) was found: 

1 
t a 

Actual Mariner deck (t^) and bottom shell (tg) thickness 
are t D = 1.12 (in) and t s = 0.815 (in). If these figures are 
taken for the c r y = 30 x lO^ materials and the other materials 
are ratioed accordingly, two similar tables of frame spacings 
may be generated from Table 4, one based on t D and one based 
on tg . 

Table 5 R atioed Plating Thicknesses 
cr y (psi) t D (in) t s (in) 

30 x 10 3 1.12 0.815 

70 x 10 3 0.48 0,35 

110 x 10 3 0.305 0.222 

150 x 10 3 0.224 0.163 
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Table 6 

Frame Spacing Based on to (In ) 

cr y /E 10 x 10 6 20 x 10 6 30 x 10 6 * 40 x 10 6 

30 x 10 3 14.85 21.00 25.80 29.60 

70 x 10 3 J*. 19 5.91 7.26 8.39 

110 X 10 3 2.12 3.00 3.66 4.25 

150 x 10 3 1.32 1.88 2.30 2.65 

Table 7 
Frame Spacing Based on t p (in ) 

30 x 10 3 20.50 29.00 35.50 40.90 

70 x 10 3 5.73 8.11 9.98 11.48 

110 x 10 3 2.91 4.11 5.01 5.81 

150 x 10 3 1.82 2.59 3.15 3.65 

The Mariner, as constructed of medium steel, used a 
frame spacing of 30 (in), about half way between the two 
frame spacings shown in the Tables above, based on t s and t D . 

Since the tables above yielded such small values of 
frame spacing for the high yield strength materials, it was 
considered not feasible to use these directly. The points 
were plotted, however? and the shape of the curves was used, 
along with the knowledge that the optimum frame spacing for 
32/30 steels in other studies was about 30 inches, to obtain 
the prediction below: (The predictions are written above the 
slash, ) 
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Table 8 
Estimate s of Fr ame S pacings (In ) 
(from curve shape of Tables 6 and 7) 

0y /E 10 x 10 6 20 x 10 6 30 x 10 6 *K) x 10 6 

30 x 10 3 17.5/I8 25/25 30/30 35/33 

70 x 10 3 10/9 18/18 23/24 28/27 

110 x 10 3 8.5/9 16/15 22/21 26/25 

150 x 10 3 7.5/9 15A5' 21/21 25/25 

In the interest of running as many runs as possible 
without changing frame spacings, these numbers were later 
rounded up or down to consolidate as much as possible. The 
numbers shown under the slash are those run as the center 
frame spacing, with one spacing three inches above and 
below the center spacing, making up a total of 48 trial runs 
in the first data set. 

In retrospect it probably would have been better to 
start at 30" and drop 5" at a time, for better consolidation 
of frame space groupings. Second and third data sets became 
necessary, since the frame spacings below 15" were not parti- 
cularly useful. In attempting to predict some starting 
point, possibly more data runs were made than would have been 
necessary if 30" were tried for all materials and then 5" 
changes were made for all materials. Figure 11 shows that 
what were finally defined as the optimum frame spacings vary 
only from 19-28 inches, 
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Definition of Optimum Frame Spacing 

A word should be said about this "optimum" frame spacing 
•that has been used so glibly. When the plots were first pre- 
pared and the curves seen for the first time, it became appa- 
rent that no definition could be found to define a true 
optimum frame spacing. 

One feels inherently, or at least the author does, that 
the optimum must be the least weight solution that is also 
the least cost solution. This feeling has no rational explan- 
ation to back it up however. One ship owner with a specific 
problem may demand the least weight solution at a sacrifice 
In cost premium. Another owner may desire the least cost 
solution and sacrifice some weight (cargo) capacity. 

Another factor complicates this optimum determination. 
At the time that the E and oy values are being varied over 
the range of the matrix all the other "variables" of these 
hypothetical materials are being held constant. (Except that 
allowable stresses vary with cr y ) This means that the density 
and the cost especially are being held constant at medium 
steel values. The decreases in weight per foot then become 
only a reflection of the reduction in the size of scantlings. 
And the cost ratios within a c v family are almost constant 
for each frame spacing. See Figures 7t 8, 9 and 10; plots of 
the hypothetical materials results in families, by yield 
strength. Since the cost of materials has been set constant 
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at the cost for medium steel, the cost ratio changes also 
reflect mainly the scantling reductions. 

When one considers the density and material costs held 
constant as above, it appears that the plot of weight ratio 
versus cost ratio really reflects only weight. The fact of 
material costs "being held constant seems to effectively de- 
couple the cost-weight relationship, This suggests that pos- 
sibly the optimum is the least weight solution. 

Never-the-less the author decided to locate the "opti- 
mum, " by definition for the purposes of this paper, at that 
point where the slope of the curve equals ^5 degrees. This 
has some sense of indicating that the least cost solution 
is of the same relative importance, to the ship owner, as 
the least weight solution. Table 10 in Appendix 1, shows 
the "optimum" frame spacing for each of the hypothetical 
materials. Figure 11 is a plot of these optimums - yield 
strength versus frame spacing, in families of constant 
modulus of elasticity. 

Presentation of Results 

Figures 12 and 13 are plots of common frame spacings 
of 21 and 24 inches, respectively. Figure Ik is composed of 
the combination of all hypothetical materials results, 
(Figures 7, 8,9 and 10) with Figure 12 superimposed. 
Another set of prediction points, similar to those of 
Figure 2 are also superimposed on Figure 1^. 
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Figure l^f- also shows a steel performance norm curve. 
This curve is the result of the studies of Reference [>l6] 
and later studies and was constructed by connecting tangents 
to Weight Ratio vs. Cost Ratio curves for steel vessels. 
These curves reflected real materials with real costs and 
real fabrication ratio estimates. As such, their results 
forced the talling-off , of the norm curve, downward and to 
the right, as yield strengths increased. This reflects the 
decreasing weights due to the higher yield strengths and 
also the increasing costs of such high strength materials. 
The results of the hypothetical materials do not show this 
"tailing-of f " tendency because all material costs and fabri- 
cation ratios for all materials was kept constant. Every 
indication for real materials, to date, points toward 
having to •'pay" a premium in cost for increased yield 
strength. This would cause "real" results to conform more 
closely the performance norm. 



79 



CHAPTER VI 
DISCUSSION AND FUTURE RECOMMENDATIONS 

General Comments 

The results of this thesis have been discussed in the 
various chapters where the detailed explanation of the 
various investigations was presented. This section will not 
recap those discussions. Instead this section will discuss 
three main areas: Mistakes made, Computer program documen- 
tation, and Recommended future studies. 

The author believes that lack of discussion of what went 
wrong is a major stort-coming of most theses. Another short- 
coming of most computer oriented theses is the program docu- 
mentation. 

M istakes Hade Along the Way 

Most theses expound, some at great length, on the suc- 
cesses achieved. Very few mention the failures that occurred 
along the way. The author believes this to be a great waste 
of valuable lessons learned, 

Many mistakes were made in handling the computer runs 
for this thesis. First and foremost, any handling of a com- 
puter program, consisting of 75 subroutines and over 6000 
Fortran cards, with data cards included, is bound to result 
in a few errors just from size alone; to say nothing of at- 
tempting to modify such a program. Several runs were complete 
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busts because of personnel handling errors at the computation 
facility. Several runs were encountered,- where more than one 
error was "reported" by the compiler and not all were found 
and corrected by the author - thus two error runs resulted 
instead of the one. 

This program is too LARGE to fit the WATFOR compiler - 
learned the hard way, but money was refunded. 

' Two limitations have been found in the present version 
of the program. These are mentioned as recommended future 
studies, also. 

Reference [16] used version 5 of this program to design 
both longitudinally framed and transversely framed, as well 
as combination or mixed framed ships. Version 6 was found 
to enter an unexplained infinite loop whenever any longitu- 
dinally framed section was attempted. 

A second looping instance occurs when the frame spacing 
is increased beyond 3.0 feet. This has occurred at 3.25 and 
3.5 feet for real material properties and for hypothetical 
materials. 

It is strongly recommended that the program be run from 
a compiled binary object deck. (this amounts to a reduction 
to about 3000 cards) Even better would be storage on data 
cells, disk or tape depending on the computer system being 
used. 
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' Computer Program Documentations - Not What They Should 
Be 

As mentioned in Chapter III, the author believ.es that 
what presently passes for program documentation is next to 
useless. The useles'sness of the diagrams and flow charts, 
that usually are called documentation, can be forcefully 
shown by examining the output of the AUTO PLOW routine now 
available. This computer routine takes any program or sub- 
routine source deck and diagrams it - or "documents it", if 
you prefer. Such "documentation" is pretty but little more 
useful than the Fortran listing itself. 

The author recommends that meaningful documentation, 
useful to the user who follows the programmer, be made a 
requirement for all computer programs. Appendix 2 is an 
example of what the author believes is the documentation 
needed; absolutely needed, by the user who wishes to modify 
a program. The longer the program the greater is this need. 

The user who must modify a program, (and this occurs 
how often at a University?) , needs to know what each variable 
is - not what kind of a diagram is used to show a DO LOOP, 
Re needs to know what are the outputs of each subroutine and 
what are the calling/called interrelations between the various 
program and subroutine levels. Most of all any user needs 
an accurate, detailed and fully descriptive set of directions 
for preparation of input data. 
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It should be acknowledged that most of the work of 
the directions for data input for this program was done ' 
previously and only modified here. 

Recommendations for Future S tudies 

The first recommendation is rather obvious after reading 
the Mistakes Made section. Investigate and fix the looping 
problems concerned with longitudinal framing and with frame 
spacings over 3.0 feet. 

Another hypothetical materials study should be conducted 
using constant, but realistic, costs for each group of equal 
high yield strength. This should cause the plots of weight 
ratio versus cost ratio to "tail-off", as the steel perform- 
ance norm curve does. It would also show much more explicit- 
ly the effect of changing the modulus of elasticity, since 
presumably some curves would fall above, and some below, 
the performance curve. 

Further studies are also needed to perfect the predic- 
tion techniqii.es attempted here. 
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APPENDIX 1 

This appendix contains tabulated 
results of the various studies conducted. 
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APPENDIX 2 

This appendix contains descriptions 
of subroutines and their interrelations and 
definitions of the program variables. 
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.DEFINITIONS OF VARIABLES 

IDIM1 = Number of decks -l- 1 (for innerbottom, if one exists) 
+ number of longitudinal bulkheads •:* 1 (for bottom 
shell) -J- number of side sections [IDIM1 = KSECT(7)] 

IDIM2 = Number of decks + 1 (for innerbottom, if one exists) 
+ 1 (for bottom shell) -J- number of longitudinal 
bulkheads [JDIM2 = KSECT(^) + KSECT(3)] 

KDIM1 - 500 divided by IDIM1 

KDIM2 = Aj-oo divided by IDIM2 

IDB11, IDIM2 S KDIM1, and KDIH2 are used to dimension arrays 

12=2, I3~3i 1^=^: used as counters or indicators 

KSHIP a Counter to indicate ship number being processed. 

KSHIP1 = Counter to Indicate ship number about to be 
processed 

BASV/T a Bs.se weight t used to form weight ratio for summary 
of results. Ratio is current ship results divided 
by base. 

BASCOS = Base cost, used to form cost ratio, as above, 

KEY ~ Print indicator - KEY=0, print all intermediate 

results and summary j KEX-1 , print only final plate 
thicknesses and summary of results, 

NM ~ Number of material sections or horizontally divided 
areas that the ships depth is divided into e 

nif = km + i 

MATRL = Material type number for all transverse structure 

TMARAT ~ Ratio of allowable a? stress for medium steel to 
allowable cro stress for material being used for 
transverse structure . 

EPR1E = Ratio of modulus of elasticity of material used for 
transverse structure to modulus of elasticity of 
steel e 

PPPRI - Ratio of quantity (1-A~) for steel to quantity 

(1-X^) for material used for transverse structure, 
where X==Poisson ratio* 
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OFRA.T - Ratio of number of open floors to number of closed 
floors per hold. 

HLRAT = Ratio of hatch length to hold length, 

BP - Length between perpinduclars [=SHIP(l)] 

B = Breadth [=SHIP(2)] 

H = Draft [=SHIP(4)] 

SALARY ~ Labor wage rate (dollars per man-hour) 

NUMUSD = Number of materials read In as data to fill 

properties of materials array (Maximum number is 5) 

NUMCOS = Number of sets of cost data read in as data. This 
number must match NUMUSD above * 

IKEY = Indicator for printing plate thicknesses; not print- 
ed on initial cycle, printed on final cycle. 

TOTCOS = Total cost per foot 

NSPACE - Number of longitudinals in a given width 

DIST = A given width of plate 

STDWTH - Standard width of plate; set as a DATA statement 
in subroutine PLTSMS , 

WIDTH = Spacing of longitudinals WIDTH = nUIcE 

ZTOP - Height of top of longitudinal bulkhead 

ZBOT = Height of bottom of longitudinal bulkhead 

ZMOD - Combined plate and stiffener section modulus 

COM BAR = Combined area, plate and stiffener 

COMBIN = Combined inertia, plate and stiffener 

P = Proportionality constant 

E = Effective width of a plate 

T = Plate thickness 
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CXHD = A.B.S. modification to head 

SY = Y coordinate 

SZ = Z coordinate 

IP = Section indicator index 

KP = Plate index 

MBR = Used to call material type number 

S3BAR = A particular tertiary stress 

S1BAR * A particular primary stress 

HEAD = A particular head 

TLAT = Thickness of plate required by lateral load 

INDEX - A particular index 

Kl = Section indicator index 

K2 = Plate index 

Z = Z coordinate 

LP = Linear program result indicator 

TOLD = Plate thickness, called old or previous 

TNEW = Output of AVERG » average of TOLD and new T in 
cycling,, 

KUP = Upper interface index 

KLO =3 Lower interface index 

SUP = Stress at upper interface 

SLO - Stress at lower interface 

SACT = A primary stress 

A = Transverse spacing 

B - Spacing between scans 

MAT - Material type number indicator 
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TBUK = Thickness required due to buckling 

BRKRAT = Bracket ratio set at 0.95 

BRAT = Calculated bracket ratio 

OK = Indicator for bracket ratio within tolerance 

ZABS = A particular section modulus 

GIRVAL = A particular section modulus 

DMIN - A.B.S. constant 

TCLAB - Total cost of labor 

TFAC = Construction time factor 
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SHIP 



KSECT 



(1) 
(2) 

(3) 
(*0 
(5) 
(6) 

(?) 

(8) 

(9) 
(10) 
(11) 

(1) 
(3) 

(2) 

(5) 
(6) 

(7) 



DEFINITIONS OF ARRAY VARIABLES 



Length 

Breadth 

Depth 

Draft 

Bilge radius 

Center vertical keel (CVK) height 

Dead rise 

Camber 

Transverse spacing 

Longitudinal spacing 

Web spacing 



No. of decks 

No. of longitudinal 

If No Inner Bottom 

(1) 

(2)+l=(l)-KL 

2*<1) 

(5)+l=2*(l)+l 

(5)+(3)=2*(i)+(3) 



bulkheads per |- section 
If an Inner Bottom 

(2) = (D+l 

W = (2)+l=(l)+2 

(5) = C0+(D»2*(l)+2 

(6) « (5)+l=2«(l)+2 

(7) - (5)+(3)=2*(i)+2*(3) 



BHD 



(1(1) =s Section No, bounding top of bulkhead I 
(2,1) =3 Section No e bounding bottom of bulkhead I 
(3 i I) - Distance of bulkhead off center line 

Longitudinal bulkheads are intercostal to 
decks. BIID(1,I) and BHD(2,I) are integers 



SPCNOM(I) 



Zero for transverse framing 
Nominal longitudinal spacing 
(I'th) section 



•i* modification for 



ZINT 



ZIP 



(1,1) = Height to lowest interface of section J 
(2,1) =s Height to highest interface of section I 

(including camber on main deck) 
(3 i I) = Material type No. 

(1,1) = Height of I'th interface 

(2,1) » Material type No, below interface I 

(3,1) = Material type No, above interface I 



NTYP(I) Type No, for deck beam support condition 

YSUP(I) Location of side girder (Coordinate off center line) 

HDSTR(I,J) 

Head load acting at stiffener location 



109 



HDSEAM(I r J) 

Head load acting at plate seam location 



PROPTX 



(1,0) - Limiting primary stress (psi) 
(2(,J) ss Limiting secondary stress (psi) 
(3»J) = Limiting tertiary stress (psi) 
(k, J) a Yield stress (psi) 
(5»J) = .Ultimate stress (psi) 
(6, J) ~ Elastic modulus (psi) 
(7»J) = Poisson ratio 
(8, J) = Specific weight (lb. /in 5 ) 
(9 i J) = Price (per pound) of plates 
(10 f J) - Price (per pound) of stiff eners 
(11, J) = Labor factor (relative to mild st€ 
J = Material type No e 



el) 



AM I TOT 



CONST 



(1) 
(2) 

(3) 

(4) 

(5) 
(6) 

(1) 
(2) 
(3) 

W 
(5) 



Area 

Moment 

Inertia 

Neutral 

Section 

Section 



axis height feet 
modulus to keel 
modulus to deck 



above keel 



Deadrise angle (radians) 

Height to bilge top 

Length from center line to 

turn 

Length of bilge turn 

Angular distance (radians) 

top along shell 



bottom of bilge 



from CVK to bilge 



GIRTH(I) Girths (lengths) of sections 



AMI 



STR 



SEAM 



1,J,K) - Area 

2,J,K) = Moment 

3,J,K) = Inertia 

1,J,K) = Stiff ener Y coordinate 

2,J,K) = Stiff ener Z coordinate 

3,J,K) = Proportionality constant 

4,J f K) = Material type No., filled by. ASGSTR 

1,J,K) - Girth (length) of seam 

2,J,K) ~ Y coordinate of seam 

3 f J,K) = Z coordinate of seam 



SPACE(I) Standard width i breaks up a width into standard 
spacing s 

KSTR(I) Number of stiff eners per section, for each section 
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HFHAT(I) 

DKHT(I) 

HDABS(I) 

VLOAD(I) 

TBEAH 



HBEAM 



BOTF 



CLFLR 



REVFR 



WEB 



WT(I) 

INI 



| of full hatch width of deck (I) 

Height of clock (I) above "base line 

Head (A.B.S.) for transverse structure for deck (I) 

Side column load (A.B.S.) or shear for deck (I) 

=; Span 

= Girder value or section modulus 

= Area 

= Side frame 

= Deck frame or outboard deck frame 

« Inbos.rd deck frame, if present 

= Half span 

- Girder value 

~ Area 

= Hatch beam or outboard hatch beam 

= Inboard hatch beam, if present 

1) = Span 

2) = Girder value 

3) ~ Area 

1) ~ Span 

2) = Thickness of CVK 

3) = Area 

1) - Span 

2) = Girder value 

3) ~ Area 

1) = Web depth 

2) - Web thickness 

3) = Flange width 
k) ~ Flange thickness 

Weight of transverse Structure per linear foot for 
section (I) 



1.*. 


i * ) 


2,«, 


"*•"* \ 


J t 1 


*>c \ 


M, 


*>$• \ 


*.2, 


-*£■ \ 


*.3i 


"VC 1 


1,*, 


JC \ 


2,* 


%;* \ 


3,* 


%j- \ 


*.i. 


"*** j 


*,2, 


•*;- ) 



(1) = 1 

(2) « 1 

(3) - 1 

(4) = 2 

(5) - 2 

(6) - 2 

(7) - b 

(8) « 5 



This array controls entry point 
after data groups are changed. 
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ERR(I) Error check array to chock all input data groups 
present on first run 

CON (1) r= /^HBASI/ 

(2) - AHDETA/ 

(3) = Ahpraj-i/ 
{k) « Ahprop/ 

(5) = Ahmate/ 

(6) = Ahstru/ 

(7) = Ahtran/ 

(8) = Ahcost/ 

(9) = Ahrun / 
(10) = Ahstop/ 

CLAB(I) Cost of labor for section (I) 

CMAT (1) = Cost of plates 

(2) ~ Cost of stiff eners 

(3) = Cost of transverse structure 
(/.j.) = Total cost of material 

CWTT(I) Converted weight of transverse structure for 
section (I) 

AP(I) Volume of plates in section (I) 

AS(I) Volume of stiff eners in section (I) 

APC(I) Modified volume, for fabrication ratioing, in 
section (I) 

ASC(I) Modified volume, for fabrication radioing, in 
section (I) 

WP(I) Weight of plates in section (I) 

WS(I) Weight of stiff eners in section (I) 

SF(I) Areas of plates in section (I) 

SCANT Standard stiffener properties 

(1) = Web depth 

(2) * Web thickness 

(3) = Flange width 

(k) ~ Flange thickness 

(5) - Neutral axis (to keel) 

(6) = Area 

(7) - Inertia 

(8) - Flange modulus 
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SIG 

NPLT(I) 
PLT 

WTLMAT 
STSPAR 



(1) = ctjl - Primary stress 

(2) - 02 - Secondary stress 

(3) ~ 03 - Tertiary stress 

Number of plates per section 



(2.*,*) = 



PLT thickness 
Material type No 



(1^ J) = Material type No, [Decimal] 
(2, J) ss Total weight for material J 
J s= Material type No e 

(1»J) =3 Bending moment (foot tons) 

(2, J) ~ Limiting stress (psi) 

(3? J) = Required section modulus - top 

(^• f J) = Required section modulus - bottom 

J » Material type No. 

The program will fill this array if data is not provided, 
However* the limiting stress will be always the same for the 
same ship dimensions and will not depend, on yield strength 
or other properties. These values are not used within the 
program but ore displayed for comparison with program 
generated results, 
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APPENDIX 3 

This appendix contains detailed 
information concerning make -up of data 
for the modified midship section program, 
version six. 
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PREPARATION P.F _DATA_INPUT 

Tne following pages present detailed instructions for 
preparation of data for computer input 6 

Each computer run must contain, as a minimum, one data 

set for each of the following groupings: 

$MS BASIC GEOMETRY 

$MS DETAIL GEOMETRY 

IMS FRAMING SYSTEMS 

|MS PROPERTIES OF MATERIALS 

$MS STRUCTURAL DATA 

IMS TRANSVERSE STRUCTURE PARAMETERS 

$MS MATERIAL ASSIGNMENT 
COST DATA 



The order of the sections is not important. Each is 
independent of the other; except that $MS COST DATA should 
always follow |MS PROPERTIES OF MATERIALS and each material 
in section $MS PROPERTIES OF MATERIALS must have a respec- 
tive cost in the section $MS COST DATA, Even if all costs 
are the same , they must be repeated to equal the number of 
materials submitted. (Maximum number of materials is 5) 

The program will test to see if all eight (8) groupings 
are filled before starting execution. If all are not filled 
a fatal error message will be generated, and execution will 
terminate. 

Once the first execution is completed, any section may 
be changed and the execution repeated with the changed data. 
No recheck or recount of the eight groupings is made, The 
program continues to execute , on the command $MS RUN, 
utilizing whatever data Is in the eight groupings? whether 
the data is the original or has been changed. 
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There are four levels of re-entry into the execution 
cycle, once the first execution has been completed, Table 16 
which follows p contains a sample data deck setup and the 
values of the control arrays that the data generates. If 
any of the following data groupings are changed the program 
must recompute the ship's geometry 1 BASIC GEOMETRY, DETAIL 
GEOMETRY, and FRAMING SYSTEMS. This is the "most expensive" 
way to re-enter the execution cycle s 

The function of the CON array is to determine what data 
group Is being read. When any of the above three data groups 
arc read the INI array is equal to one and entry point 
indicator (Nl) is set equal to one* (Ml = 1) The entry 
point indicator is always set equal to the lowest value of 
the INI array that is encountered* On the first run this is 
naturally set to one e 

The next level of entry is used when any of the follow- 
ing data groups are changed? PROPERTIES OP MATERIALS, 
MATERIAL ASSIGNMENT, and STRUCTURAL DATA. In this case the 
geometry does not have to be recalculated. These are the 
wore useful data groups to change for parametric studies and 
they are "cheaper" to change than the first three, When 
TRANSVERSE STRUCTURE PARAMETERS or' COST DATA are changed, 
even less re-calculation is necessary. These are the "cheap- 
est" changes to iTiake c 

This "relative cost of change" of the various groupings 
should be borne in mind when setting up the sequencing of 
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data runs. If frame spacings arc to be varied, for instance, 
the number of changes should be Minimized. This means that 
all runs with common frame spacings should be run consecu- 
tively} changing material assignments, material properties 
or cost data, etc. 

The following pages detail the identification of sec- 
tions and the input formats for the various data groupings. 
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IDENTIFICATION OF SEC TIONS 

Number all decks starting -with main deck, number inner 
bottom next after last deck, number the bottom shell next 
(this section follows the bottom shell shape up to the center 
vertical keel height or inner bottom height) , number the side 
shell from inner bottom to lowest deck next, number the re- 
maining side shells between decks next, number any longitu- 
dinal bulkheads next (from lowest to highest) 



Section 1 
Section 2 
Section 3 
Section 4 
Section 5 
Section 6 
Section 7 
Section 8 
Section 9 



Deck 1 (main deck) 

Deck 2 

Deck 3 

Inner Bottom 

Bottom Shell 

Side Shell, inner bottom to next upper deck 

Side Shell, next higher side shell 

Side Shell, next higher side shell 

Longitudinal bulkhead 



Refer to Figure 15, which follows, for illustration, 
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$MS BASIC GEOMETRY 



Guide to preparing basic geometry data input 
(See subroutine BASIC) 



Card 1 



Columns 1-10 
Columns 1.1-20 
Columns ?.1_30 
Columns 'Jl-io 



Length ( ft ) 
Breadth (ft) 
Depth (ft) 
Draft (ft) 



[Format will be shown here; the words integer 
or decimal will be used along with the stand- 
ard Fortran notation: Decimal, ^FlO.ty] 

Card 2 

Column 5; Number of decks (less than 9, do not count 
inner bottom) 
[Integer, 15] 

Card 3 

Inner bottom information: 

a) If an inner bottom is to be included write; 
INNER BOTTOM, beginning in column 7 

"b) If no inner bottom is to be included, write; 
NO INNER BOTTOM, beginning in column 7 

Card ^ 

Column 5; Number of longitudinal bulkheads in half 
section 
[Integer, 15] 



Card 5 



Columns l-io 

Columns 11-20 

Columns 21 -30 

Columns 31-tyo 



Bilge radius (ft) 

Center vertical keel (CVK) height (ft) 

Deadrise (ft) 

Cambsr (ft) 



[Decimal, ^FIO.2] 

Card 6 

Data for forming weight and cost ratios and print 

indicator KEY, 

Columns 1-10; Base weight (tons/ft) 

Columns 11-20; Base cost ((?/ft) 

Column 22; Print KEY-- (one or zero) 

KEY=0 t Print all detail; KEY=1 , Print summary 

only, 

[Decimal 2F10.2, Integer 12] 
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GEOMETRY 



Guide to preparing decks hatch and longitudinal 
bulkhead detail geometry (See subroutine DETAIL) 

One card for each deck, in same order as numbered. 
Columns 1-10; Deck height above baseline (ft) 
Columns 11-20; Total hatch width, if any (ft) 
[Decimal 2P10.2] 

One card for each longitudinal bulkhead; longitudinal bulk- 
heads are intercostal to decks 

Column 5 5 Section number bounding top of bulkhead 
Column 10; Section number bounding bottom of bulk- 
head 
Columns 11-20; Distance of bulkhead off centerline (ft) 
[Integer 215, Decimal F10.2] 
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PB AMINO SYSTEMS 



Guide to preparing framing system data input 
(See subroutine FRAME) 



Card 1 

Columns 1 
Columns 11 
Columns 21 



-10 $ Transverse frame spacing (ft) 

-20? Nominal longitudinal spacing (ft) 

-30; Web spacing; if ship is longitudinally 

framed (ft) 
[Decimal 3^10 „ 2] 



One card for 
Column 



Columns 8 



each section 

7; Framing system, type: 
Longitudinally framed 
Transversely framed 
Modification of nominal spacing 
or «0,6 j ft added or subtracted 
nominal spacing) or leave blank 
modification. 
[3X, Alpha A4 f Decimal P13.2] 



L; 

T; 
=20; 



(ie: + 
from 
for no 



0.3 
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plS PROPERTIES OF MATERIALS 

Guide to preparing properties of materials data 
(See subroutine PROM AT) 

Card 1 

Column 7; Number of materials to be read in. 
[Integer I?] 

One card 2 and one card. 3 for each material 

Card 2 

Column 7; Material type number 
Columns 8-17; Yield strength (psi) 
Columns 18-27 j Ultimate strength (psi) 
Columns 28-37 ; Allowable stress factor 
[Integer 17, Decimal 3F10.0] 

Allowable stress factor, type: 

0,0; program generates average of allowable 

stresses ratioed from medium steel 

values, based on yield and ultimate 

strengths. 
Some other factor (positive); program will 

multiply medium steel values by factor. 

Medium steel values are: o-, =20000 (psi) 

02=27000 (psi) 



Card 3 



Columns 1-12 
Columns 13-22 
Columns 23-32 



Oo=32000 (psi) 



Modulus of elasticity (psi) 

Poisson ratio 

Specific weight (lbs/in^) 



[Decimal F12.0, 2F10.0] 
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STRUCTURAL DATA 



Guide to preparing structural data input 
(See subroutine STRUCT) 



Card 1 

Column 

Columns 



7; Number of different materials for 
which input data is provided. 
8-17; Bending moment (ft-tons) 
[Integer 17 s Decimal F10.2] 



One card for each material numbered above 

Columns 1-10? Limit on highest primary stress allow- 
able (psi) 
Columns 11-20; Require^ section modulus at deck 

(in 2 -ft 2 ) 
Columns 21-30 j Required section modulus at bottom 

(in 2 -ft 2 ) 
[Decimal 3F10.2] 

Note: 1, This program will fill all the above data 

slots if none are submitted , You must enter 
1 in column 7 of the first card and put a 
blank card in for card two* If the program 
computes the above data, it is based on 
medium steel formulas * and on ship dimensions 
Therefore all material sections would give 
the same results, and only one will be 
printed, 

2 e This data is not used internally but is 
displayed for comparison with computer 
generated results* 

3, For materials other than steel the minimum 
recommended data input is: a) Fill card l t 
b) Fill limit stress (columns 1-10) on a 
card for each material r c) Program will fill 
remainder, which are simple division and 
unit changes, 

k. If you desire to submit one required section 
modulus, for both top deck and bottom . place 
it in columns 11-20. 
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TRANSVERSE STRUCTURE PARAMETERS 



Guide to preparing transverse structure parameters 
input data (See subroutine TRANSV) 



Card 1 

Columns 1-10 j Ratio of open to closed floors per hold 

(ies 20 open and 5 closed-ratio-^. 0) 
Columns 11-20; Ratio of hatch length to hold length 

(le: 60 ft. hold, 30 ft. hatch-ratlo=0.5) 
[Decimal 2F10.2] 



One 



card for 
Columns 



each deck 
5; Type 



of support condition for deck thru- 
beams (See following illustrations) 
Columns 6-2^ j Location of side support girder, if one 

exists (given as distance off centerline 8 
in feet) (If none, leave blank) 
[Interger 15 f Decimal F19.2] 
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$M£L .-MATER I AL.. AS S IGME_NT 



Guide to preparing material assignment data input 
(See subroutine HATASG) 

Card 1 

Column 7; Number of material areas that ship's depth 

is divided into. 
Column l'+; Material type number for all transverse 
structure, 
[Integer 21?] 

One card for each material area, numbered above. 

Column 7; Material type number for area 
Columns 8-20; Height to lower interface of area (ft) 
Columns 21-30; Height to upper interface of area (ft) 
[Integer 17, Decimal F13.2, F10.2] 

Note: l e The lowest interface must be 0.0 at the 
keel for the first area. 

2. The uppermost interface height, at the 
main deck, must include camber, if any. 

3. The upper interface of the previous area 
always becomes the lower interface of 
the next area above , 

4. Cards must be ordered, lowest area to 
highest, keel to main deck. 
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$MS COST DATA 



Guide to preparing cost input data 
(See subroutine INCOST) 

Card 1 

Column 7; Number of materials for which cost data is 
provided 
[Integer 17] 

Note: This number must be the same as the number 
of materials for which data is provided in 
data group £MS PROPERTIES OF MATERIALS 
(ie: each material must have an individual 
set of cost data) 

One card for each material numbered above 

Column 7; Material type number 
Columns 8-17; Cost of plates (Dollars per pound) 
Columns 18-27; Cost of stiff eners (Dollars per pound) 
Columns 28-37; Fabrication factor - Based on costs of 

fabrication of medium steel (ie: 1,2=20^ 

greater fabrication cost) 
[Integer 17, Decimal 3F10.3] 

Last card 

Columns 1-10; Labor wage rate (Dollars per man-hour) 
[Decimal F10.2] 



130 



APPENDIX k 

This appendix contains a copy 
of a sample input data dec!:. 
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APPENDIX 5 
This appendix contains sample 
output of the computer program. The 
results are output generated by the 
input data shown in Appendix k. The 
results of the first and last "ship" 
are included. 
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SHIP 1 SUMMARY 



MATERIAL SECTION ASSIGNEE MS 

TOP , 3.0 SICE , 3.0 BOTTOM , 2.C 



LCNGITUDIN^L WEIGHTS 
MATERIAL TYPE NO. 
3. 



WEIGHT ( TC.NS PER FCCT) 



2 .2 0' 



TRANSVERSE WEIbHT 

MATERIAL TYPE NO. 

2,0 



WEIGHT ( TCNS PER FOOT ) 



TOTAL WEIGHT 



749 



COSTS (DOLLARS PER F03T) 

COST OF MATERIALS = 629.9 

CCST CF LABCR = 2116,6 

TOTAL COST = 2746.5 



WEIGHT RATIO = 0.601 COST RATIO = 1.C23 
RATIOS BASED ON 



BASE WEIGHT - 6.236 

BASE COST - 2 6 64.7 
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SECTION 



INERTIA 



D E T A I L S 



NEURAL AXI S HEIGHT 
SECTICf^ NCCULUS 



19*50 FEET ABCVE KEEL 



TO KEEL = 21010.79 SC.IN.-FT. 

TO DECK = 15455.44 SG.IN.-FT. 
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PLATE 



THICKNESSES 



NC. 



THICKNESS 



FET/L TYPE 



CECK PLTS 





(PLT. 


SECTION 


1 




1 






0. 


41 


2 






0. 


4 1 


3 






0. 


42 


4 






0. 


42 



CECK PL1S 





(PLT . 


SECT ION 


2 


> 


1 






0. 


32 


2 






0. 


32 


3 






0. 


32 


4 






0. 


32 



DECK PLT5. 



CFLT. StCTICN 



0. 


3 


0. 


30 


c. 


3C 


0, 


30 



1^9 



NC 



THICKNESS 



METAL TYPE 



lf\f\ER BOTTOM PLTS. 



(PLT. SECTION 



) 



1 
2 
3 
4 
5 
6 
7 
8 



0.32 
0.31 
0.31 
0.31 
0. 31 
0.31 
0.31 
C.32 



?. 
3 
3 
3 
3 

■3 

3 
3 



SI- ELL PLTS . 



1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 



(PLT. SECTION 



) 



0.38 
C.35 
0.35 
0. 35 
0.3 5 
0.3 4 
0. 31 
0.3C 
0.31 
C.31 
0.31 
0.31 
0.31 
0.34 
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SHIP 



S U P M /i R Y 



MATERIAL SECTICiN ASSIGNMENTS 

TOP , 3.0 SICE , 2.0 BOTTOM , 3.C 



LONG IT ID IN AL WEIGHTS 
MATERIAL TYPE NO. 

2.0 

3. 

TRANSVERSE WEIGHT 

MATER UL TYPE HO. 

1.0 

TOTAL WEIGHT 



WEIGHT (TCNS PER FOOT) 
1. 447 
2.54 9 



WEIGHT (TCNS PER FOOT) 
1.266 

5.262 



COSTS (uCLLARS PER FCCT) 

CCST CF MATERIALS = 1516.7 
COST OF LABOR = 1867.5 

TOTAL' COST = 3386.3 

WEIGHT RATIO = 0.844 CCST PATIO = 1.261 



RATIOS 8 AS ED CN 



BASE WEIGHT = 6. 2^6 
BASE CCST = 2634.7 



* v- v $ # *>- :;- * * y: 



^^^^^^^-^^^^^^^^.^.^^^^^ 



'- ^ * V 7 =r } * i v ^ A i) i; 



15? 



L i T £ 



T ri I C 



\ c S S E S 



NC. 



THKfcNESi 



1 - L 1 > 9 



CECK FLTS. 



( 3 LT. SECTION 



1 ) 



0.61 
C. 6? 
0.62 
0. ~3 



CECK PLTS. 



(FLT. ScCTIJT, 



1 

2 

- 



0.55 
0, 55 
C. 55 



CECK FLTS. 



(PLT. ScCTIC\ 



3 ) 



0,49 
C. 47 
0.4 7 
0.47 
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NO. 



THICKNESS 



FETAL TYPE 



INNER BOTTCM PLTS. 

(FLT. SECTION 4 ) 



1 
2 
3 

4 
5 
6 
7 

8 



C.57 
0.56 
0.56 
C.56 
0.56 
C.5b 
.56 
0.57 



SHELL PLTS. 



(PLT. 


SECTION 5 


1 


C.66 


2 


0.53 


3 


0. 5c 


4 


0.5 8 


5 


0.5d 


6 


0.53 


7 


C. 5 7 


8 


0.56 


5 


0.66 


10 


0.66 


11 


.7 1 


12 


C. 71 


13 


C.85 


14 


0.91 
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SECTION 



INERTIA 



DETAILS 



NEURAL AXIS HEIGHT 
SECTICN fCCULUS 



19.67 FEET £BCVE KEEL 



• • • • 



TC KEEL 
TO DECK 



36495.64 SO, IN. -FT. 
27275.13 SG.IN.~FT. 
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